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Abstract 

We summarize and extend previous results on the comparison of threshold resummation, performed, using soft- 
collinear effective theory (SCET), in the Becher-Neubert approach, to the standard perturbative QCD formalism based 
on factorization and resummation of Mellin moments of partonic cross sections. We show that the logarithmic accu- 
racy of this SCET result can be extended by half a logarithmic order, thereby bringing it in full agreement with the 
standard QCD result if a suitable choice is made for the soft scale fi s which characterizes the SCET result. We provide 
a master formula relating the two approaches for other scale choices. We then show that with the Becher-Neubert scale 
choice the Landau pole, which in the perturbative QCD approach is usually removed through power- or exponentially 
suppressed terms, in the SCET approach is removed by logarithmically subleading terms which break factorization. 
Such terms may become leading for generic choices of parton distributions, and are always leading when resummation 
is used far enough from the hadronic threshold. 

Keywords: QCD, soft-gluon, threshold, resummation, soft-collinear effective theory 



1. Soft resummation and scale choices 

Threshold resummation [ 1 , 2 1 plays an important role 
in extending and stabilizing the accuracy of perturba- 
tive results, and it may be of some relevance even for 
hadronic processes which are quite far from thresh- 
old |3] 01, due to the fact that the underlying par- 
tonic process can be rather closer to threshold than 
the hadronic one 0. All-order resummed results are 
known to lead to a divergent series when expanded out 
in powers of the strong coupling: this is physically due 
to the fact that resummation is obtained by choosing as 
a scale of the parton-level process the maximum energy 
of the radiated partons E |7), which tends to zero in 
the threshold limit. The divergence can be tamed by in- 
troducing suitable subleading contributions, such as ex- 
ponentially suppressed terms outside the physical kine- 
matic region [8|, or power-suppressed terms l9l [l0ll . 



'Speaker. 



In Ref. ifTTl it was suggested, within the context of a 
SCET approach to threshold resummation, that the di- 
vergence can be tamed by making a hadronic choice of 
scale. In SCET this is possible because resummed re- 
sults are characterized by a "soft scale" fi s : the Becher- 
Neubert (BN) scale choice consists of expressing fi s in 
terms of kinematic variables of the hadronic scattering 
process. The meaning of this choice is not obvious in 
the conventional QCD approach, where, because of per- 
turbative factorization, the partonic cross section, which 
is being resummed, is independent of the hadronic kine- 
matic variables. 

In Ref. lfT3l we have clarified this issue by explic- 
itly exhibiting a relation between the ji/ s dependent re- 
summed SCET result, and the standard (fj. s independent) 
QCD expression. Specializing to the BN scale choice 
(while taking the Drell-Yan process lfT2l as an exam- 
ple) we were able to show that in the SCET result, with 
the BN scale choice, the divergence is removed through 
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QCD: 


A(a s ) 


D(a s ) 


go(a s ) 


accuracy: a" In* N 


SCET: 


Ecusp(^j) 


Jw(a s ) 


H, iDY 


accuracy: a" \n k (jj s /M) 


LL 


1-loop 




tree-level 


k — 2n 


NLL* 


2-loop 


1-loop 


tree-level 


2n - 1 < k < 2n 


NLL 


2-loop 


1 -loop 


1 -loop 


2n-2<k<2n 


NNLL* 


3 -loop 


2-loop 


1 -loop 


2n - 3 < k < 2n 


NNLL 


3 -loop 


2-loop 


2-loop 


2n-4<k<2n 



Table 1: Orders of logarithmic approximations and accuracy of the predicted logarithms in perturbative QCD (first header) and SCET (second 
header). The last columns refers to the content of the coefficient function. In Ref. [12], only the N*LL* counting is considered for the SCET 
resummation. 



terms which are logarithmically subleading in compar- 
ison to the logarithmic accuracy of the SCET result of 
Ref. Ifl2l . which is by half a logarithmic order lower 
than that of the standard QCD result. 

Here we show that the accuracy of the SCET result 
of Ref. Ifl2ll can actually be increased to the same level 
as that of the QCD result, and we rederive, within this 
higher accuracy, our master formula, which thus be- 
comes particularly transparent. We then use this im- 
proved master formula to discuss various problems re- 
lated to the BN scale choice. 



2. Resummation of the Drell-Yan process in SCET 
and QCD 

We consider for definiteness inclusive Drell-Yan pro- 
duction, but the same discussion applies to other pro- 
cesses, such as Higgs production in gluon-gluon fusion, 
with minimal modifications. The dimensionless invari- 
ant mass distribution <x(t, M 2 ) = — ^Sf, with M the 

v ' 1 rtr dM- ' 

invariant mass of the pair and cr the leading order par- 
tonic cross section, can be written schematically (omit- 
ting a sum over partons) in factorized form as 



ct(t,M 2 ) = J *i;Qc(z,M 2 ), 



(1) 



where _£ is the parton luminosity, s the hadronic center- 
of-mass energy squared, and t = — , so that the 
hadronic threshold limit is r — > 1. The perturbatively 
computable coefficient function C(z, M 2 ) is normalized 
so that C(z, M 2 ) = (5(1 - z) at leading order in the strong 
coupling a s . In the sequel, without significant loss of 
generality, we shall always choose the renormalization 
and factorization scales jj. 2 F = [i 2 R = M 2 . 

Standard QCD resummation follows from a Mellin- 
space renormalization-group argument J6] |7) : indeed, 
at the resummed level, both the convolution Eq. ([1]) and 



the gluon radiation phase space factorize, so at the re- 
summed level one may write 

cr(N, M 2 ) = f cIt t n - 1 cr(r, M 2 ) = £(N) C(N, M 2 ) 
Jo 

(2) 

where the A^-space resummed coefficient function has 
the form [DEI 

C QCD (JV, M 2 ) = go (a s (M 2 )) exp S {m 2 , ^ J (3) 



where 



SIM 2 



M 2 
N 2 



dz — 

o 1 -z 



(4) 



Jm 1 



1?- 



2A(a s (M 2 )) + D(a s ([l - z] 2 M 2 )) 



The functions go(a s ), A(a s ) and D(a s ) are power series 
in a s , with £ (0) = 1 and A(0) = D(0) = 0. 

Because resummation is obtained through exponen- 
tiation, it might seem natural to also exponentiate the 
function go. However, unlike S(M 2 , M 2 IN 2 ), go is inde- 
pendent of and only depends on a s (M 2 ). As a conse- 
quence, it turns out that simply including an extra term 
in £0 at each order increases the logarithmic accuracy of 
the coefficient function by half a logarithmic order. This 
is summarized in Table [T[ where the logarithmic accu- 
racy obtained by including a given number of terms in 
^o, A, and D is summarized. A given accuracy means 
that all and only the logarithmically enhanced contribu- 
tions to the coefficient function listed in the last column 
are correctly predicted. At leading logarithmic (LL) ac- 
curacy only the largest power of In at each order in a s 
is predicted; adding one order in each of the functions 
A, D, and go one then obtains the next-to-leading log- 
arithmic (NLL) accuracy, which correctly predicts two 
powers more, and so on to N^LL accuracy. However, if 
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£0 is exponentiated and a power counting is performed 
at the level of exponents, it may seem more natural to 
include one less order in £0- This results in the N A LL* 
accuracy, also shown in table, which is lower by one 
power of InN at each order in a s than the N A LL accu- 
racy. 

The resummed SCET expression for Drell-Yan pair 
production is given by lfT2lP ] 



C SCET (z,M 2 ,^) = //(M z )t/(M z ,^)5(z,M 2 ,^) (5) 

where H(M 2 ) (hard function) is a power series in 
a s (M 2 ), 

S{z,M 2 ,£) = Sdy In — + F ,J \ — - 

(6) 



(soft function) depends on 

"A J„2 



1 = f ~T" r ™^p (a s (M 2 )) . r cusp (o' J ) = A(a s 



). 

(7) 

soY(L,fi s ) is a series in a s (fi 2 ) with L-dependent coeffi- 
cients, and 



In U(M 2 ,i4) 



t 

Jm 2 



dfi 2 



- cusp 



(a,0u 2 ))ln 



M 2 



Jw 



(8) 



where yw(a s ) is also a series in a s , with jw(0) = 0. The 
scale fi s is a soft matching scale of the effective theory, 
and C SCET formally does not depend on it, up to sublead- 
ing terms. However, the SCET result resums powers of 
In 'g, so this choice of scale determines what is being 
resummed. 

Again, a given logarithmic accuracy is obtained by 
including a finite number of terms in the perturbative ex- 
pansion of the functions which determine the resummed 
result, namely r cusp , yw, H and Sdy, according to Ta- 
ble [l] In Ref. Ifl2ll . only the N*LL* accuracy was con- 
sidered: in fact, the order called N A LL* in Tab. [T] is 
actually referred to N A LL Ref. fl2l . which might be 
the source of some confusion. Computations using ei- 
ther of these two definitions of the logarithmic accuracy 
have been presented in the past, either in the contex of 
QCD (see e.g. Ref fl4l . where N A LL* is referred to as 
N*LLi n fl) or SCET (see e.g. Ref. El, where N A LL is 
referred to as N A LL')- 



'The resummed expression as given in Ref. 1 12) actually depends 
on several hard energy scales, which here for simplicity are all taken 
to be equal to the hard scale M 2 . 



Here we show that in fact the N A LL* of Ref. Ifl2ll can 
be promoted to higher N A LL accuracy, by inclusion of 
the terms listed in the table. This result is obtained in 
the next Section, by explicitly computing the relation 
between this improved version of the SCET result, and 
the QCD result. 

3. Comparison at NNLL 

An analytic comparison between the QCD and SCET 
resummation formalisms can be performed [13] in 
space, where the QCD result is naturally constructed, 
and where it admits a convergent perturbative expan- 
sion in powers of a s . Namely, we determine the ratio 
C r (N, M 2 ,/4) between the QCD and SCET expressions, 
Eqs. (|3) and @: 

C QCD (N,M 2 ) = C r (N,M 2 ,iJ 2 )C SCET (N,M 2 ,n 2 ). (9) 

In Ref. im we computed C r {N,M 2 ,^i 2 s ) to NNLL, us- 
ing the definition of NNLL of Ref. 03, which in Tab.[TJ 
we call NNLL*. Here we show that the accuracy of the 
SCET expression can be upgraded, and the comparison 
can be carried out at full NNLL according to the defini- 
tion of Table [TJ 

We first rewrite the QCD result in the more conve- 
nient form irm 



C QCD W M 2 ) = g (a s (M 2 )) exp^ QCD (M 
with N = Ne y and 
lo(a s ) = 8 o(a s ) exp 



#2 



(10) 



2( 2 A(a s )+l^j3 —a 2 

J 7T 



(ID 



M--S)-Jr¥h^: 



M 1 



A(a s ) = —a s + —a s + — a* +0(a s ), 



, (12) 
(13) 



£>(a s )= -D(a s )-2b—/3 a 2 s 
2 n 



D 2 a 2 + 0(a]), (14) 



and/Jo = (HQ - 2n f )/(12n). 

The functions A, D and £0 are computed at NNLL 
order according to Table [TJ the relevant coefficients can 
be found in Ref. lfT3l . except the two-loop contribution 
to £0, which can be determined by matching the expan- 
sion of Eq. ([3} to the NNLO Drell-Yan cross section in 
Ref. M, 



fo(a s ) = l + -C F (2<r 2 -4) 

7T 



M. Bonvini et al. /Nuclear Physics B Proceedings Supplement 00 (2013) 



4 



n 2 16 



592 

<^|--<r 2 2 + ^2 + 28<r 3 - 



12 
5~ 
72 



+ Cf I y£- 70f 2 -60&+ 4 



112 127 
+ "/|86-- 9 -6 + — 



1535 \ 
^2~j 
511 



(15) 



and the coefficient A 2 , which we give here for complete- 
ness 



4C F 



67 



•2& C 



10 



(16) 



The explicit expression of A3 is not needed here. 

On the other hand, at NNLL the Mellin transform of 
the SCET result, Eq. |5]), can be written as Ifl3l 



C SCET {N,M\rf) = H(M 2 )E\—,ijl 2 s 



M 2 



with 



H(M 2 ) = H(M 2 ) exp 



x exp o SCET (N,M 2 ,n% (17) 



E [^r-H\ 2 N 



2 n 
M 2 2 \ 

eX P 



(18) 



"2 ;r 



(19) 



>Sscet(N, 



,M 2 ,rf) = 

Jm 2 



dfi 2 



M 1 i 1 



, sp (a s (n 2 )) In 



fi 2 N 2 



7w(a s ) = 7w(a s ) - % —/3 a 2 , 
2 n 



(20) 
(21) 



and fwfas) = Diets) at this order. 

In comparison to Ref. [13], we now also include 
the two-loop contributions to the functions H and Sdy, 
which were given explicitly in Ref. |[T2l . Note that, in 
order to be accurate to order a 2 , the definition of the 
function E slightly differs from Ref. iTHl . 

Putting everything together we find 



C r {N,M 2 ,ri]) 



go{a s (M 2 )) I M 2 

; — exp<S \u„ 

H(M 2 )E(§, M l) V s N 2 



(22) 



with 



-/ 2 m 2 \ r Ml ^ 2 dn 2 \ t 2 x 



M 1 



D(a s (M 2 )) 



(23) 



It is easy to see that 

go (a s (M 2 )) 



l+0{a)(M 2 )), (24) 



H(M 2 )E(M 2 , M 2 ) 
so to NNLL accuracy Eq. ( p2| i can be written 

C r (N,M 2 ,fi 2 ) 



2 2 , E(M 2 ,M 2 ) $ 2 M 2 " 



(25) 



Using the 2-loop expression of suy from Ref. Ifl2l in 
Eq. ( fT9l >, we find 

e|^,^J= l+E l (L) as ( M 2 ) + E 2 (L)al(jj 2 ) + 0(al) 

(26) 



where 



£,(L) = ^ 2 , 



(27) 



A 2 
Y2 
CaCf 



1 , L 3 Ai L 2 A 2 

£ 2 (L) = — —L -Bo + + LD 2 

2V ' 128 3 8 ^ 2 16 



607 67 



11 



324 + 144 ^ 4^ 72 ^ 3 



162 72 t2 36 



and 



L = In ■ 



M 

'1W 2 



(28) 



(29) 



Note that L = when the two arguments of E are equal 
to each other. 

Eq. ( |25| establishes our first new result. Indeed, it 
is immediate to check that C r (N,M 2 ,n 2 ) = 1 for /j s = 
M/jV, up to subleading (NNNLL*) terms. This means 
that with this scale choice the SCET result now repro- 
duces the QCD result to full NNLL accuracy, rather than 
to the lower NNLL* accuracy of Ref. lfT2l . 

We are however interested in studying C r for generic 
scale choices, and in particular with the BN scale 
choice. The result becomes especially transparent by 
casting the ratio E(M 2 , M 2 )/E {%,H 2 S ) Eq. ((B) in the 
form of an exponential of an integral, of the same kind 
as the form adopted in Eq. ( |23] l. This can be done at 
the price of including terms of order a] or higher in the 
ratio, which is allowed at NNLL. The ensuing expres- 
sion of C r is particularly simple and suitable for analytic 
comparisons. It should however be kept in mind that a 
numerical comparison of the SCET and QCD expres- 
sions should rather be performed using the exact expres- 
sion Eq. ( p5j ), and possibly also retaining the subleading 
terms in Eq. (24 1. 
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We get 

E(M 2 ,M 2 ) 



In 



+ a 2 s (ji 2 s ) 



A, L j 



A, L 2 

4~t 

AjL 2 

16 2 



D 2 L 



+ 0(a]). 



Using 



A+l n 



k 



M IN 1 ^,,2 



d/j 2 , k M L 
H 2 N 2 



(30) 



(31) 



and taking the running of a s into account, we finally 
obtain 



2 , r M /A ' dn 2 m 2 

A x a s (ii 2 ) A 2 a 2 (fi 2 ) 



X 



A(a s (// 2 )) 



16 



(32) 



which is our NNLL master QCD-SCET comparison for- 
mula. It generalizes to full NNLL the result of Ref. lfT3i 
Its most notable feature, which determines the relative 
accuracy of the comparison, is that (recall the expansion 



Eq. ( fT3"] i) the exponent in Eq. (32 1 is of order a^. Note 



that this is however due to the exponentiation Eq. ( 30 1. 
If one does not exponentiate (as in the original SCET 
expression), when expanding C, in powers of a s , terms 
proportional to A; and A 2 only cancel up to 0(a 2 ). 



where C QCD (z, M 2 ) and C SC et(z, M 2 ,jj 2 ) are the inverse 
Mellin transforms of Eqs. ( fTOj i and ( fTTj i, respectively, 
and in the SCET case after performing the inverse 
Mellin transform at fixed jj s we have set yU, = M(l - r). 
Of course, C QCD (z, M 2 ) is given by a divergent series in 
powers of a s (M 2 ), so it should be understood as the 
order-by-order Mellin inversion up to arbitrarily high 
but finite order. Using Eq. |9} we find 

<r QCD (T,M 2 ) = j *cr SCET (pM 2 ) C r (z,M 2 ,M 2 (l-T) 2 ) 

(35) 

where C r (z, M 2 , M 2 (l - t) 2 ) is the inverse Mellin trans- 
form of Eq. ( p5| l, performed at fixed fi s and evaluated at 
H S = M{\ - t). 

In order to compute C r (z, M 2 , M 2 (l - t) 2 ) it is conve- 
nient to rewrite Eq. ( |25| ) as 

, , E(M 2 ,M 2 ) £0u 2 // 2 ) -/ , M 2 \ 
C r {N,M ,£) = , , ,\\ exp<S U 2 — . 

(36) 

The first ratio is just a function of a s (M 2 ) and a s (fi 2 ), 
independent of N, while 



expSirf,^) = l+F r (a,Q&,L), (37) 



where F r (a s , L) is of order a 3 s . 

The inverse Mellin transform can be now computed 
using the results of Ref. |[T3l . We find 



4. The Becher-Neubert scale choice 

As briefly discussed in Sect. [T] the BN approach 
is based on the idea of choosing for fi s a scale de- 
termined by hadronic, rather than partonic kinematics, 
namely /j s = M(l - t). In Ref. IfTTl [T2l a more gen- 
eral choice yu., = M(l - t)#(t) was considered, with 
g(r) = const. +0(1 - t): the distinction may be relevant 
for phenomenology, but it is immaterial for our present 
goal, which is to determine the logarithmic accuracy of 
the SCET result with this scale choice. 

Since the variable t refers to hadron kinematics, the 
comparison can only be performed at the level of the 
physical cross section, Eq. ([T}. We therefore define 



cr QCD (r, M ) 



o- scet (t, M 2 ) 



Jt Z "^(z) 



C QCD (z,M 2 ), 



(33) 



C SCET (z,M 2 ,M 2 (l -r) 2 ), 
(34) 



cr 0CD (r,M 2 ) = 



E(M 2 ,M 2 ) 



EM,*® 



ct scet (t, M 2 ) 



d 



+ F, .,».(//•;). 2. ) Kr..V/\^ 



where 



£(t,M 2 ,£) ! ' 



e rfr(£) 



(38) 



(39) 



for /i s = M( 1 -t). We have shown in App. B of Ref. Ifl3l 
that£(r,M 2 ,£) can be expressed in terms of derivatives 
of <t scet (t, M 2 ) with respect to ln(l - r), up to terms 
suppressed by positive powers of (1 - r): 



E(t,M 2 ,£) = £q(£) 



k=i) 



d k o- SCET (j, M 2 ) 
d\n\l-T) 



[l+OCl-T)], 



(40) 

with coefficients which do not depend on t. It follows 
that the term proportional to F r in Eq. ( f3"8] > does not con- 
tain any extra logarithmic enhancement with respect to 
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<t scet (t, M 2 ). On the other hand 
E(M Z , M 2 ) 



= 1 + E 2 (0) [a 2 (M 2 ) - a 2 (^ 2 )] + ... (41) 
= 1 + 4E 2 (0)/3 a 3 s (M 2 ) ln(l - r) + 0(a 4 s ). 



We conclude that the leading difference between the 
QCD and SCET expressions is 

cr QCD (r,M 2 ) - cr SCET (r,M 2 ) = 

= ct scet (t, M 2 ) [a^4A)£ 2 (0) ln(l - r) + . . .] , (42) 

where the ellipse denotes terms which are either of rel- 
ative order O(a^), but without any logarithmic enhance- 
ment, or 0{a*). 

Because the log counting is now done at the level of 
hadronic cross sections, it is based on counting powers 
of ln(l -r). Also, because the SCET result violates stan- 
dard QCD factorization (i.e., it does not factorize upon 
Mellin transformation), the difference between cr QCD and 
<x SCET depends on the parton luminosity (it is not univer- 
sal) through cr SCET itself. A generic leading-log term in 
<x SCET has the form 



K*ln 2 * + "(1 -t), 



(43) 



where ln 2k (l - t) is due to the leading log behavior of 
the coefficient function, and ln p (l - t) generally comes 
from the parton luminosity. 

If we assume that the parton luminosity does not lead 
to any logarithmic enhancement, then 



^"OCD C"si 



jt+3 j n 2it+l 



(1-t) = a';in 2 ^ 5 (l-T), (44) 



where we have set h = k + 3. This corresponds to a 
NNNLL* correction. It is interesting to observe that, 



had we used the exponentiated version Eq. ( 32 1 of C r , a 
NNNLL, rather than NNNLL* correction, would have 
been obtained. The argument can be generalized to the 
case in which C, is computed to all orders in a s rather 
than just to order a 2 . Indeed, no leading logarithmic en- 
hancement arises from the factor 1 + F r (a s (/-i 2 ), L); the 
only possible source of powers of ln(l - r) in C r is the 
ratio E(M 2 , M 2 )/ E(fi 2 , fi 2 ). It is easy to see, however, 
that all terms in the expansion Eq. ( [41) are at most of 
order a k s {M 2 ) ln 4 ~ 2 (l - r), with k > 3. Thus, the conclu- 
sion Eq. ( |44| holds to all orders in a s . 

In conclusion, we restate three observations which 
were already made in Ref. 11131 . First, we note that 
the BN scale choice removes the divergence of the 
perturbative expansion at the cost of introducing loga- 
rithmically suppressed non-universal terms. This is to 



be contrasted with the commonly used Minimal pre- 
scription [8 1, which also introduces non-universal terms 
(with support outside the physically accessible kine- 
matic region) but are more suppressed than any power, 
or with the Borel prescription l9l [T0ll . which introduces 
power-suppressed but universal terms. 

Second, we observe that quite in general we do ex- 
pect PDFs to contain logarithmically enhanced terms. 
In this case the terms introduced by the BN scale choice 
to tame the perturbative divergence can become lead- 
ing or even super-leading (i.e., more logarithmically en- 
hanced than the leading log). 

Finally, we remark that threshold resummation is of- 
ten useful in situations where r is far from threshold, but 
nevertheless the partonic subprocess is close to thresh- 
old 13] ID . In this case M(l - r) ~ M, and conse- 



quently C, Eq. ( 25 1 is actually leading log. 



The phenomenological implications of our results re- 
main to be investigated. They are potentially of con- 
siderable interest, given the increasingly important role 
that threshold resummation, in its various implementa- 
tions, is playing for LHC phenomenology. 
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